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ABSTRACT 



Ion acoustic waves have been artificially excited and detected in 
the Steady State Plasma Facility at the Naval Postgraduate School by a 
magnetic excitation device and an optical detection system. The dis- 
persion of these waves has been studied as a function of the exciting 
frequency and the confining static magnetic field. 

Dispersion curves are presented which compare the experimental data 
to a new dispersion equation which is derived in the present thesis. 
After the validity of this dispersion relation is demonstrated, the 
measurements are used as a diagnostic tool, and five operating para- 
meters of the plasma facility are deduced from the measurements. 

Theories are presented to explain the behavior of the dispersion 
curves and three of the five operating parameters. 
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Section 1 



INTRODUCTION 

For the purposes of the present thesis a plasma may be defined as 
a quasi-neutral electrically conducting mixture of positive ions and 
electrons plus a mamber of neutral background gas atoms. Due to the 
particular plasma production mechanism employed in the Steady State 
Plasma Facility of the Naval Postgraduate School, the electrons are at 
a relatively high temperature (22 electron volts) on the cylindrical 
beam axis, while the ions have thermal energies of the order of .09 ev 
and the neutral gas atom energies are approximately .01 ev. 

Due to the high electrostatic restoring forces created by even a 
small internal charge imbalance, a plasma has the ability to maintain 
overall internal charge neutrality. If a charge imbalance should 
develop due to the random thermal motion of the electrons and ions, the 
particles will rearrange themselves to shield the remainder of the 
plasma from the resulting electric field. Thus beyond a certain average 
distance from the imbalance the electric field so created will not be 
felt- This distance is called the Debye length and is, of course, a 
function of the plasma parameters. In a true plasma the Debye length 
must be much less than the smallest dimension of the plasma system. 

The Debye length is givei 




where 



€o is the permittivity of free space 



K is the Boltzmann constant 



f is the electron temperature 
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N is the equilibrium plasma density 



e is the electron charge 



Ion acoustic waves are longitudinal plasma density oscillations 
occurring at frequencies well below those of the electron plasma oscil- 
lations. The plasma ions and electrons move exactly in phase only at 



the frequencies employed in the present investigation, however, the ap- 
proximation is an excellent one. The ion oscillations are slightly 
larger in amplitude than those of the electrons. The ratio of ion per- 
turbation velocity to electron perturbation velocity can be shown to be 



Due to their small mass and high mobility the electrons see these 
low frequency waves as quasi-static perturbations. Ion acoustic waves 
are similar to ordinary sound waves in a neutral gas, except that the 
propagation mechanism is the Coulomb interaction between charged particles 
rather than the short range collision forces between molecules. Their 
existence was first predicted by Tonks and Langmuir in 1929 based on the 
fluid analysis of a two component plasma, 

Before proceeding to an account of previous work, the Tonks- 
Langmuir dispersion relation for ion acoustic waves will be derived. 

The following assumptions are made: 



Lt) = 0; for ^ ^ 0 in-phase motion is only approximately realized. For 



2 

shown rnar ^ ) 

2Tf 




where Ml is the wavenumber. In Section 2 it will be 
10 , so it may be neglected compared to 1. 



1) The plasma electrons have negligible momen^.tum compared to 



the plasma ions. 



2) The plasma ions have zero temperature and are singly 



charged. Thus there is no ion pressure. 



3) There is always 100% ionization within the plasma. 
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4) There is no electric field associated with the wave. 

5) The Debye length is much smaller than the wavelength, 

that is Xp« • 

6) The drift velocity of the ions is zero. 

7) The unperturbed plasma has uniform density. 

Under these assumptions the momentum transfer equation for the ions can 
be written as: 

+ p) = + V •'^][ (N + n)v] - (N + n)VXz 



where 



M is the ion mass 

p is the electron pressure perturbation 
P is the equilibrium electron pressure 
V is the ion velocity inducted by the wave 
N is the equilibrium plasma density of charged particles 
n is the plasma density perturbation of charged particles 
e is the electron charge 
Bq is the static magnetic field strength 
z is a unit vector parallel to the direction of 
Subtracting the momentum transfer equation for the ions in the un- 
perturbed system from the above equation and linearizing to the limit of 
small perturbations gives: 



^ '^p = (MV) -IlNV0r +JlNVr.e 



M 



where is the ion cyclotron frequency 



( 1 ) 



eBr 



M 



Now for perturbations of the form ^ equation (1) can 

be solved to give: 
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and 



NVr. = 



NVq = 



( 

M Tr ^ 



il2 _o2 



-g / 1 

M Tr“ \g_2 _<o 2 



Therefore 



V ■ (BVe? - NV,S)= ^ . 4- I 4 , ) 



Taking the divergence of equation (1) gives: 



or 



^^^p = (NV) - S2^- (NV0? - NVi-e) , 



(NV) + — ^)(- 

M ' M ^j.2 r d r 



.2 -4^2 



-) (2) 



But by the linearized equation of continuity, 

'^n 



• (NV) 



yr 



= iCvin 



Therefore equation (2) becomes: 



^2p = _ O ^ (.V£_ ^ J_^)(_ 



M "2>r2 ' ^2 _^2 

Assuming there is no radial electron temperature variation, 



and 



2/"<i2n 1 ^ _.\2, 



c^C 









JTL ^n ^ 1 *^n N 

^ r2 r 5r 



where 



2 - 

M 






(def init ion) 



and 7* is the electron temperature. 



12 



Rearranging and defining 




this equation becomes: 

x^n* * + xn’ + x^n = 0 

Its solution is; 

n - n^e 

where Jq is the zero order Bessel function. 

Therefore it must be true that 

= (2.405.. .)2 

where the boundary condition that the perturbations vanish at the 
position r = ^ is used. (The use of this type of boundary condition 
will be more fully discussed in Section 2.) This can be rewritten as; 



^ (\2,SL^ 



5.79 









Ai2 



- 1 ) 



Now for an infinite medium an 



^2 



or 



i.= 






(3) 



Thus for an infinite medium the application of a static magnetic field 
has no effect. Equation (3) is the well known Tonks -Langmuir dispersion 
relation for ion acoustic waves. 

From an elementary acoustic analogy the ”hot” electrons provide the 
perturbing and restoring pressure gradients, while the massive ions pro- 
vide the inertia. This dispersion relation predicts a phase velocity of 

= c = ( ^ )^. Also, since , the group velocity equals 

the phase velocity and, there is no dispersion. 

A similar dispersion relation will be derived in greater detail 
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including additional effects in Section 2. 

Ion acoustic waves can be self-excited or externally excited both 
in gaseous discharges of low percentage ionization and in highly ionized 
plasmas. The existence of self-excited waves in low density gaseous 
discharges has been observed by Alexeff and Neidigh |^4j and by Little and 
Jones 1^5 j in the form of standing waves in the discharge tubes employed. 
The present investigation deals exclusively with propagating waves. Two 
earlier experiments (one on standing waves and one on propagating waves) 
will be described. 

Ion acoustic waves have been excited by electrostatic methods (grids 
and probes |^7j) and by magnetic perturbations of the local electron 
diffusion coefficient The magnetic perturbation method was used 

in this experiment. Its theory of operation appears in Section 3. 

Once the ion acoustic waves have been excited, they can be detected 
by grids, probes, and light sensitive detectors. In this experiment both 
probe detection and photomultiplier detection were used, although only 
the photomultiplier detection system was used to make quantitative mea- 
surements. Detector theory and operation are also covered in Section 3. 

Wavelength and attenuation measurements were made as functions of 
exciting frequency for various static magnetic field strengths. The 
validity of the dispersion relation is demonstrated, the data are anal- 
yzed and the results presented as dispersion and attenuation curves. 

The new dispersion and attenuation theories derived in Section 2 
contain 5 adjustable parameters. Comparison of the experimental data 
with theory allows computation of these 5 parameters, whose values are 
presented as functions of the static magnetic field strength. The ion 
acoustic waves are thus employed as a diagnostic tool. Physical inter- 
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pretations of the experimental results (dispersion curves, attenuation 
curves, and parameter behavior) are also given in Section 5. 



Wong, Motley, and D'Angelo have carried out experiments on the dis- 
persion of ion acoustic waves in cesium and potassium plasmas in strong 
magnetic fields, 

Fried and Gould have shown that the Vlasov equation predicts 

ion acoustic wave propagation even in the absence of collisions in the 
plasma. Under this condition, however, the analysis predicts that the 
waves will be strongly damped due to a "trapping" of ions moving with 
velocities close to the wave phase velocity. This can occur if there are 
ions whose temperatures are close to the electron temperature (T = r>- 
From equation (3) the wave phase velocity is c = ( )^, and if T = 7" 

there will be ions meeting this condition. 



This type of collisionless damping was first predicted by Landau. 

It was observed by Wong, Motley, and D'Angelo, and their work provided 
the first experimental proof of its existence. The conditions of their 
experiment were such that ion-neutral collisions were negligible (the 
percentage ionization varied between 40% and 907o) . Their plasma was 
produced by surface ionization of cesium and potassium atoms on a hot 
tungsten plate. This method guaranteed that electron and ion tempera- 
tures were approximately equal. The damping was so strong that standing 
waves were not observed. 

The Q-3 Machine (Princeton) was used for the experiment. The fea- 
tures essential to the experiment appear in Fig. 1. The plasma was pro- 
duced at one end of the machine, allowed to drift longitudinally, and 
terminated approximately 90 cm away by a similar tungsten plate. Radial 
confinement was achieved by adjusting the longitudinal magnetic field to 
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FIG. I 
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12,000 gauss for cesium and 6000 gauss for potassium (these values were 
chosen to avoid possible ion cyclotron resonance). 

In the Princeton experiment an electrostatic excitation method was 
used, in constrast to the magnetic excitation method employed in the pre- 
sent investigation. The perturbing signal was applied to the excitation 
grid and modulated the percentage ion transmission through the grid. The 
movable receiving grid detected the density perturbations. 

Phase velocity calculations were made by measuring the slope of the 
curve of phase delay vs. receiving grid position (the movable grid) for 
the propagating wave. One result of the experiment was that the phase 
velocity was found to be independent of frequency for the range from 15 
kc to 100 kc (ion cyclotron frequency was 140 kc for cesium and 240 kc 
for potassium). This conforms to the predictions of the dispersion re- 
lation derived in the present thesis, although the frequency range rela- 
tive to the ion cyclotron frequency is different. 

The damping observed by Wong, Motley, and D'Angelo was found to be 
in excellent agreement with the theory of Landau damping. In Section 2 
the criteria necessary for the presence of Landau damping will be dis- 
cussed, and it will be shown that it can be ignored in the present inves- 
tigat ion. 

Figs. 2 and 3 show phase velocities for the ion acoustic waves in 
the Q-3 experiment for propagation parallel and anti-parallel to the plas- 
ma drift velocity. 

The drift velocity of the plasma was estimated from the difference 
between the phase velocities of the waves for parallel and anti-parallel 
propagation. In the present investigation the plasma drift velocity is 
one of the 5 adjustable parameters. The magnitude of the drift velocity 
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PLASMAS AS A FUNCTION OF FREQUENCY. THE DIFFERENT SYMBOLS REFER 
TO DATA TAKEN AT DIFFERENT ION DENSITIES: X« l.8xlO'° cm-3} 
t*6.0xlo'®cm"3i o*l5xlo'°cm-3i +«6.0xl0*® cm'S; 30x10*® cm-3. 



2x10 




K 
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FREQUENCY (kc/sec) 

mi 

PHASE VELOCITY OF THE UPSTREAM WAVES IN CESIUM AND POTASSIUM PLASMAS AS A 
FUNCITION OF FREQUENCY. PLASMA DENSITY IS 2xlO" cm”*. 



is machine calculated from the dispersion data. Both parallel and anti- 
parallel propagation at 1000 gauss were studied in the present investi- 
gation, and the difference in phase velocities for the two directions of 
propagation provided an independent determination of the drift velocity. 
The method is the same as that employed by W^ng, Motley, and D'Angelo. 

This determination is discussed in Section 5. 

Little and Jones have made dispersion measurements on externally 
excited ion acoustic waves in the positive column of a low pressure mer- 
cury arc. The percentage ionization in the discharge was approximately 

1%, the current through the mercury vapor discharge was fixed at 9 amps, 
and the axial magnetic field strength was 45 gauss. Their excitation 
mechanism was the same as that employed in the present investigation. 
Essentially the magnetic excitation device employed created sufficient 
variations in the local electron diffusion coefficient to modulate the 
plasma density. (This method of excitation will be explained in detail 
in Section 3.) 

Little and Jones detected the waves by observing variations in the 
light output of the plasma ^due to the plasma density variations) by 
means of photomultiplier detection. This method was also used in the 
present investigation and also will be explained in Section 3. Fig. 4 
s*hows a schematic diagram of the detection apparatus used by Little and 
Jones . 

Their experimental dispersion curves show a dependence of phase 
velocity on frequency and a low frequency cutoff below which propagation 
cannot occur. These results are not directly comparable to those of 
Wong, MDtley, and D'Angelo, since the frequency ranges relative to the 
ion cyclotron frequency are entirely different. Fig. 5 shows a plot of 
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stationary 
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FIG. 4 

DETECTION CIRCUIT SCHEMATIC 
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versus for 45 gauss external field. The frequency cutoff condition 
can be seen. 

• The results reported by Little and Jones conform to the theory de- 
rived in Section 2, although their agreement with their theory is not as 
good as in the present thesis. This is probably due to the unclear 
nature of the boundary conditions in the discharge tube of their experi- 
ment (as they have pointed out). The boundary conditions applicable to 
the present investigation will be discussed in detail in Section 2, and 
it will be shown that they can be unambiguously specified. 

The present investigation examines the dispersion of ion acoustic 
waves in a high density plasma of high degree of ionization (approxi- 
mately 90%) and considers the combined effects of ion drift, damping, 
ion temperature, and a confining magnetic field together with finite beam 
geometry. 

The present thesis has three purposes. First, an improved disper- 
sion relation is derived for a cylindrical plasma including the afore- 
mentioned effects Second, this dispersion relation is tested against 
experiment. Third, from the comparison of theory with experiment, the 
above 5 parameters for the Steady State Plasma Facility at the Naval Post- 
graduate School are computed by least squares computer search to give 
the best fit to the data. The validity of the deduced values of the 
parameters is also demonstrated. 

The derived theoretical solution to the wave equation is a complex 
equation. The real part of the solution predicts the dispersion rela- 
tion, that is wave number as a function of frequency. The imaginary part 
predicts the damping as a function of frequency. 

The solution to the wave equation is derived by assuming the ion and 
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electron pressures are scalars. This allows truncation of the set of 
moments of the Boltzmann equation after the first 2 moments, which are 
the equation of continuity and the momentum transfer equation respective- 
ly. These equations are combined with the plasma equation of state. As- 
sumption of the azimuthal and longitudinal forms of the wave behavior 
allows the resulting partial differential equation to be reduced to an 
ordinary differential equation involving only the radial behavior. This 
is solved in series form and yields the solution to the dispersion prob- 
lem when the boundary condition is inserted. 

Experimentally the waves are excited by perturbing the local elect- 
ron diffusion coefficient perpendicular to the confining static magnetic 
field. Increasing the local electron cross-field diffusion coefficient 
decreased the local electron density. Since the ions and electrons move 
in phase, this creates a deficiency in the local plasma density of ions 
and electrons. This density perturbation propagates parallel to the 
static magnetic field. In the present investigation the frequency of 
the applied perturbations was varied from 187 kc to 352 kc . 

The wavelength and amplitude of the ion acoustic waves are measured 
as functions of frequency. The results are reduced by machine computa- 
tion and fitted against theory. As has been previously mentioned the 5 
adjustable parameters in the theory are ion drift velocity, the effective 
ion collision frequency for momentum loss, ion thermal energy, plasma 
beam diameter, and the e-folding distance for the radial density profile. 
The first 3 are characteristics of the plasma facility, while the last 2 
are functions of the applied static magnetic field. Independent deter- 
minations of 4 of these 5 quantities have been made, and these values are 
compared to the results of the least squares adjusted values obtained in 
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the present investigation. The results are presented as functions of the 
static magnetic field, and theories are presented to explain their be- 
havior. 
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Section 2 



BASIC THEORY AND DERIVATION OF THE DISPERSION RELATION 

If Landau damping is to be a factor in wave propagation, the dis- 
persion analysis must begin with the consideration of the particle dis- 
tribution functions, since the mechanism involved is the resonant or 
near-resonant transfer of energy between the wave and those particles 
having thermal speeds approximating the wave phase velocity. If there 
are more particles whose thermal speeds are slightly less the wave phase 
velocity than there are particles whose thermal speeds are slightly 
greater than the wave phase velocity, there will be more particles taking 
energy from the wave than there are giving up energy to it. Then there 
will be a net loss of wave energy, and Landau damping will exist. In 
any system having a Maxwellian distribution there will be some particles 
having thermal speeds approximating the wave phase velocity, but for the 
particular plasma under consideration this is a negligible fraction of 
the total number of particles. 

As will be shown later, the phase velocity of these waves is greater 

than or equal to Spectroscopic measurements by R. Booth and 

M 

R. L. Kelly of this laboratory have established that ion energies are 
less than 0.5 ev when argon is used in the plasma facility, j^ioj Probe 
measurements made by D. M. Gall have established that electron energies 
along the cylindrical beam axis are 22 ev. ^llj Thus it may be safely 
assumed that the electron temperature greatly exceeds the ion tempera- 
ture. (This assumption is also consistent with the results of the data 
reduction presented in Section 4.) Thus there will be an insignificant 
number of positive ions capable of causing Landau damping. 

An estimate of the Landau damping may be made by examining a dis- 
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persion relation derived by Fried for the case T»T. ^12 j 






m 




1- ■ 




- 


41t 2 




r 0(1 + J 










’] J 



where 



m is the electron mass 



and 



M is the ion mass, 



For the plasma studied the Debye length is of the order of 0.02 mm, 

>-14 

-nr- i-'- 

4 

and may be neglected compared to 1. Making this approximation, 



^2 ^^2 ^ 

whereas typical wavelengths varied from 5 cm to 20 cm. Thus ^ >5^ 10 






Thus it can be seen that an upper limit on the Landau damping term (the 
imaginary term in brackets) is ( = 2.3 x 10""^. This is sufficient- 
ly small that it may be ignored compared to the damping arising from ion- 
neutral collisions. This neglect of Landau damping is even more fully 
justified if ion-ion collisions (ion fluid viscosity) are considered. 

The effect would be to introduce disruptive collisions which would des- 
troy the phase relationship between the resonant ions and the wave. This 
critical phase relationship could also be destroyed by ion-neutral colli- 
sions, but not appreciably by ion-electron collisions due to the negli- 
gible electron-to-ion momentum transfer. It should also be noted that 
the thermal speeds of the electrons are much higher than the wave phase 
velocities encountered in the present investigation. Thus the electrons 
see these waves as quasi-static and do not contribute any Landau damp- 
ing. The process of Landau damping therefore need not be considered for 
the purposes of this thesis. The derivation will neglect moments of the 
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Boltzmann equation higher than the second. 

For the conditions in this experiment under which ion acoustic waves 
propagate, ions and electrons move in phase due to the strong electro- 
static restoring forces, so the instantaneous ion and electron densities 
are equal. Although the Coulomb interaction is the mechanism by which 
these waves propagate, there is no macroscopic electric field associated 
with the waves. This assumption is justified at any frequency if the 



magnetic pressure, (* 



Br 






-) , is much greater than the plasma pressure 

B^2 



r 1 B ^ 

(essentially NKT). [13] In this experiment ( — 2 — 4000 joules/m^, 
while NAc'Y = 10 joules/m^. The following derivation will therefore omit 
consideration of any electric field associated with the waves. 

In the following derivation the quantities W, Bq, Y , T are assumed 
to be functions neither of position nor of time. 

W is the equilibrium ion drift velocity 

is the effective ion collision frequency for momentum loss 
due to ion-neutral collisions 

These assumptions are justified in the cases of Bq and T . They are 
probably approximately true in the case of T as well. Since W is defined 
as an equilibrium drift velocity, it must be independent of time, and for 
simplicity it is assumed to be independent of position as well. 

The momentum transfer equation for ions is first linearized and then 
combined with the linearized equation of continuity and the plasma equa- 
tion of state. The usual form is assumed for the azimuthal and longitu- 
dinal forms of the wave solutions, and a series solution of the resulting 
ordinary differential equation of the ion density perturbations is ob- 
tained. Insertion of the radial boundary condition yields the disper- 
sion relation which is then compared to the experimental data. 
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Assuming the electron and ion pressures are scalar quantities, the 
momentum transfer equation for the ions under the influence of the per- 
turbation is: 

■^^(P +TT + p + 'll) = + (V + W) (N + n)(V + W)j (4) 



where 



+ (N + n) (V + W) 7 + n) (V + W)X^ 



P is the equilibrium electron pressure 
T is the equilibrium ion pressure 
p is the electron pressure perturbation 
Tf is the ion pressure perturbation 
V is the ion velocity induced by the wave 
n is the ion density perturbation 
and the positive z direction is taken in the direction of Bq. 

The momentum transfer equation for the ions in the unperturbed 
plasma column is: 

^ 7 (P +ir ) = [ J [ Nw] + NW •?- _Q NWXz (5) 



Expanding equations (4) and (5), subtracting (5) from (4), and 
assuming the density perturbation is small compared to the unperturbed 
density allows the resulting equation to be partially linearized to the 
first order in small quantities. 



7(P + T ) = N + v 

+ V • y (N + n) (V + W) 



( 6 ) 



+ W • (NV + nV + nW) 
+ ^ (NV + nW) - II NVXz 
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Since solutions of the form sought, the follow- 

ing operators may be replaced by the listed variables; 






-iiO ( cd assumed real) 

iA (A complex longitudinal wavenumber) 



^ iyU complex azimuthal wavenumber) 

Also since W is assumed to be in the z direction, W = Wz and 



w -V 



w 






The radial form of the density perturbations is left unspecified for the 
present . 

Linearizing the expression V • V (N + n) (V + W) to the first order 
in small quantities gives: 

V • ^ (N + n) (V + W) = W nV0 + i'iffiV^) , 



where 



\ ^ ^ 
V = V^r + VqG + V^z. 



Thus (6) becomes: 



-1 . _ . _ A 

~j:j'7(p +10 = -i«-JNV - i«^nW + WVj. z 



+ i>^W(NV + nW) +7 n(V + W) - SINVqV +J1 NVj.9 



(7) 



From (7) are immediately obtained: 



■^(p +V)(i^- Uw -.;)) + -#:?■) 



MNVg = 



MNV^ = 



(idJ- -9 )2 +_^2 

(i^- iiw -9) - —i^JlCp +1t) 

(ii^- ijgw -t>)2 +1^ 



(7a) 



(7b) 
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(7c) 



iJ^ 



NV = 
z 



„ (p +'/!) - + iiw^n +v?WN + WV_-^^ 

M ^ c? r 

( - *) ) 



Forming the vector difference of equations (7a) and (7b) and taking 
the divergence of the resulting quantity gives: 



V . (NVer - NV^e) = ^ 



(163- iJ&J - p)2 +^2 



( 8 ) 



Neglecting volume ionization and recombination the equation of 
continuity for the ions is: 

• (N + n)(V + W) = - ’ (N + n) = itOn (9) 

Taking the divergence of equation (7) gives: 



-1 

M 



■v^(p 



+ W) = -iO 






NV + V • nW| 






+ 




NV + 




NV + V . 




( 10 ) 




NVgr - NVr.0 I 



Linearizing equation (9) gives: 

V . (N + n) (V + W) = . NV + "7 • nW (11) 



Using equation (11) to evaluate V • (NV) + • (NW) 

\ s 

and using equation (8) to evaluate y . (NVgr - NV^-O), 
equation (10) becomes: 
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( 12 ) 



1 (p +')!') = - cS^n + cJ^Wn - icOpn 
M 






■g2 

M 



(9^^ - y;7-)(p -^r) 

(i<0- 



w 



MN 



^ r 



' [ iAUc^- iiw - 9)-^ 






r 



■]<p 



+ '7f) 



(i^O - ij8.W - p +_^'- 



Under the assumption of the spatial invariance of Y T equation 
(12) becomes: 

/ 2 /-,2^ fl ^ ir (i^- iJ^ - 7 a2 ? 

lA i^W - p) + 



(13) 



.2 n2^ W 

" ^ ^ — TF 



(iO- ijg.w - ?)^ + SL 



where c 



2 , 



+ cd^n -6ijSwn - i^Pn = 0 

(definition) 

M 



and C" 



KT 

M 



(def init ion) . 



The differential equation (13) would be directly reducible to 
Bessel’s equation of order ^ if were equal to zero (that is if the 

unperturbed column had a uniform radial density distribution). This, 
however, is not the case in the present investigation. 

The mean free path for the ion-neutral collision process is greater 

_ o 

than 1 m, while the ion gyro radius is less than 10 m. Up to several 
ion gyro radii away from the beam axis the analysis of Pfirsch and 
Biermann predicts that for a Maxwellian velocity distribution and a 
Gaussian distribution of the guiding centers of the charged particles 
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there will be a Gaussian density profile. ^14j Previous probe analysis 

by Gall and Oleson has fitted a Gaussian radial density profile to N, 

at least to a good approximation, ^isj Thus N = ^ , so that 

-1— and equation (13) becomes: 

N c) r ^ 



(c" + C^) 



2^ fr ^ ^n ^ 1 S n 



m 



r Sr 






(ic^- i^W - 
(iO- ilW -9)^ 






- (c2 + c2)(2Wr,j2), 



i>6r^^(i<^ - -■?) + 

(iti- i4,W 



(14) 



+^^n -oAwh - i«vi^?n = 0 
Defining s e \ r 



and 



V 2 _ r -‘•IF + ii? ) 1 



+ 2^S^J2 

(tO-J^W + i9)^ 



equation (14) becomes 
^s2 ^ s 



s2 . s 



I!_[ s - 2sY-^V J + n/s2 -. 2 ) = 

s L (cO-AW + ip ) I / 



(15) 



Assuming a solution of the form 

n = n^.S.^ 2^ aj^j (a^ / 0) 
iVtf j=o 

and inserting this into equation (15) leads to the recursion relation 



P (j + 






"2j " TaJT 



2g 



2 '""' 2J jl (j *1^1): p (. 



|a»I 



and 



2g 



-) 






(where a^ has been chosen as — p-i ). 

2<^' l^ll 
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Thus the solution of equation (15) is: 



cc 

n(s) = n El 
j=o 



r (.i + b)s 

2 2^ j: (j + |/a| )i r (b) 



( 16 ) 



where 






b = Iris - 1 



2g 



and 

g ^ 26 ^ Aw 

(cO)-iw + iP)X2 



Equation (16) reduces to the Bessel function for the case g = 0 
(uniform density in the unperturbed column). For this case the disper- 
sion relation would assume the simple form X /^ = 2.405 where the boundary 
condition has been imposed that r = is the radial position at which the 
perturbation vanishes. The validity of this type of boundary condition 
was confirmed by the probe measurement discussed in detail in Section 3. 
Independent confirmation that this type of boundary condition is appro- 
priate can be obtained from an examination of the radial electron temper- 
ature profile as obtained by Gall and Oleson. jl5j This curve is repro- 
duced as Fig. 18 in Section 5. It is evident from an examination of 
Fig. 18 that the electron temperature drops sharply (or at least becomes 
small compared to the value at r = 0) at the position r = 7,8 mm. An 
examination of equation (14) under the condition that both c and C be- 
come small reveals that the magnitude of the density perturbation, n, 
must also approach zero. The use of this type of boundary condition is 
particularly convenient, since no assumption need be made concerning 
the existence of either a pressure node or a velocity node at the plasma 
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boundary. This problem has plagued investigators working with low ion- 
ization gaseous discharges, 

For the case g ^ 0 the dispersion relation for = 0 is: 



00 

z 



gJ P (j - -^) j 

2Jj: j: r (- -^) 



0 



(17) 



where now 

g = 28 ^A>\J 



and 



] - - 'X 


^ + iP) 


[ (6i- jg.w + i )2 J[_ 


^ c2 + c2 J 



Equation (17), or more accurately the recursion relation between 
the coefficients, is compared to the experimental data by a least 
squares direct search computer program on a CDC 1604. 
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Section 3 



EXPERIMENTAL ARRANGEMENT 

A schematic illustration of the Steady State Plasma Facility at the 
Naval Postgraduate School appears in Fig. 6. The plasma is produced by 
a hollow cathode discharge, which has been investigated by Lidsky et. al. 
j^lbj The details of its construction appear in Fig. 7. The facility is 
capable of operation on helium, neon, and argon, although only argon was 
employed in the present experiment. The neutral argon gas is supplied 
through a micrometer valve assembly to the hollow cathode. An r-f volt- 
age is supplied between the cathode and the first anode in order to ig- 
nite the arc. Once the arc has been struck, the discharge is self- 
sustaining, and the r-f voltage is turned off. 

The plasma so produced drifts through the hole in the first anode 
and into the experimental tube. This section is a 3 meter pyrex tube of 
inside diameter 10 cm. It is axially aligned with the confining magnetic 
channel to within 1 mm throughout the length of the experimental region. 

The optimum operating conditions for the facility are an arc cur- 
rent of 100 amps and a voltage of 90 volts between the cathode and first 
anode. This condition is most conveniently achieved by regulation of 
the argon gas pressure through the micrometer valve assembly. 

Differential pumping of neutral gas atoms by the diffusion pumps 
maintains a neutral gas pressure of approximately 3 x 10"^ mm Hg through- 
out the experimental region. 

The strength of the confining magnetic field is continuously variable 
from 500 gauss to 10,000 gauss. Some experimental data were taken at a 
field strength of 600 gauss, but the agreement with theory was quite poor. 
It was suspected that some drift occurred in the current supplying the 
field coils, so that the confining magnetic field was not constant with 
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SCHEMATIC DIAGRAM OF PLASMA FACILITY 
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FI6. 7 

CATHODE-ANODE ASSEMBLY AND VIEWS OF BOTH ANODES 



time (a necessary condition for taking a complete set of dispersion data). 
These data are not reported, and the minimum field strength employed was 
800 gauss. 

Upon reaching the end of the experimental region, the plasma ions 
and electrons strike the floating anode (so-called because its potential 
is allowed to float) and are recombined. They are then pumped away as 
neutrals by the vacuum pumps in the burial chamber. A photograph of this 
section appears in Fig. 8. A side view of the entire facility appears 
in Fig. 9, and a view of the cathode-first anode region appears in Fig. 

10 . 
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FIG. 9 

.SlOE VI EJ OF FLFvSM\ FACILITY 
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In order to excite the wave, a short r-f coil was tightly wrapped 
around the pyrex tube as shown in Fig. 6. It was designed for optimum 
coupling to the plasma at the center of the band of exciting frequencies 
employed. This coupling effectiveness was determined by the length of 
the coil, which was approximately that of the wavelengths excited. Ex-- 
perimentally it was found that the r-f field strength dropped appreciably 
in the last % of the coil at each end, so that these outer regions were 
ineffective in producing a perturbing field. The coil itself was approx- 
imately 15 cm long, leaving approximately 7 cm for the effective field- 
producing region. This was in the center of the band of wavelengths 



The coil was the inductive arm of a parallel '^almost resonant” tank 



to resonance at the applied transmitter frequency. This slight departure 
from the resonance condition was necessary so that the transmitter would 
see a load impedance similar to that of the radio transmission antenna 
which it was designed to feed. 

The transmitter used to power the tank circuit was a U. S. Navy 
type TAB-7 manufactured by the West ingho use Corp. It is capable of 
producing a peak r-f power level of 2 kw over the frequency range of 
100 kc to 500 kc. This full power condition was not used due to the 
inability of the r-f exciting coil to dissipate the heat generated at 
this condition. Typical operating conditions were a transmitter output 
current of 3-4 amps. The Q value of the tank circuit was estimated to 
be about 10 at the frequencies used. Thus the actual exciting current 
in the coil was approximately 30-40 amps. Direct measurement of the 
r-f magnetic fields produced yielded values of from 125 to 195 gauss. 



exc ited . 



circuit, which was tuned (with 
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depending on transmitter power. As previously mentioned, this perturbing 
magnetic field produced variations in the local electron diffusion co- 
efficient, creating plasma density variations. 

Assuming the magnetic perturbation, B, is of the form B = B^sin o)T 
and that (which was true in all cases), only the influence of the 

perturbing field on the electrons need be considered, since the massive 
ions are unable to respond directly to the rapidly varying magnetic field, 
while the electrons can instantaneously follow its variations. The dif- 
fusion rate of electrons perpendicular to the magnetic field is: 



Db 



1 + 



D 

e(B + Bp) 

m TV 



■1 



where Dg is the local electron diffusion coefficient with both B and 
Bq applied 

D is the local electron diffusion coefficient with zero applied 
magnetic field 
e is the electron charge 
m is the electron mass 

is the electron-neutral collision frequency 

eB 

Now under the conditions of the experiment, ( 1. Therefore 

m 

Dg = D \ — ° ^ 1 ^ or Dg<^ j- . 

^ L e(B + Bo) J ^ (B + Bq)^ 

Now for the general diffusion problem Dg ^ = Cy (source) 

(under conditions of steady state). Replacing by — , where 

b r2 

a is an average effective diffusion mean free path, there results 

1 o 

Thus NoC — — and (N + n) oC (B + B^)"^ or 



=Q or N = 

a ^ 



Db 



Dr 
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. o Bisin cdt 9 

(N + n) oC (Bisin<L>f+ Bq)'‘- Since B^<<Bq, (N + n) = N(1 + g 

. 2B3^sin c0^ • j. 

= N(1 + r Therefore n = 2N — — sinc-ir, and the plasma den- 

sity is modulated at the applied frequency, the strength of the perturba- 
tion being proportional to that of the r-f magnetic field. These per- 
turbations \^?ere studied over the frequency range of 187 kc to 352 kc and 
for static magnetic fields from 800 to 2200 gauss. 

Once the waves have been generated, they propagate parallel to the 
static magnetic field in both directions. Only propagation anti-parallel 
to the static magnetic field was investigated in detail. (One set of 
parallel propagation measurements was made; this allowed an independent 
determination of the ion drift velocity, W. ) 

Previous spectroscopic measurements have established that virtually 
all the light emitted by the plasma is from excited neutral atoms. |^loj 
Therefore the emitted light intensity at any point along the plasma col- 
umn is directly dependent on the local electron density. (An ion- 
electron light production mechanism would produce a light intensity xAich 
was proportional to the square of the local charged particle density.) 

Thus a photomultiplier tube can be used to detect the presence of the 
waves. The detection circuit schematic is given in Fig. 11. The photo- 
multiplier tube was housed in a light-tight electrostatic and magnetic 
shield manufactured by the Quanta Corp. A 1P21 photomultiplier tube was 
mounted inside this shield and dynode voltages were supplied by batteries. 
The entire assembly was located in a region of minimum magnetic field 
(less than 25 gauss). The optical signal was fed to the phototube 
through a 10 ft flexible unaligned optical fiber bundle manufactured by 
Optics Technology, Inc. This fiber bundle caused a 677o loss in light 
intensity, but the signal strength was large enough that this caused no 
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A Keithley Model 110 Wideband Amplifier with an overall gain of 
unity was used as an impedance matching device between the phototube 
and the Hewlett-Packard Model 450A Wideband Amplifier (which was set to 
a 20 db gain) . In order to improve discrimination against noise the 
signal was then fed into a Spencer- Kennedy Electronic Filter lyfodel 302, 
which served as a band pass filter centered at the transmitter frequency. 

The next element in the detector circuit was a Hewlett-Packard Model 
400H Vacuum Tube Voltmeter which was used to make amplitude (attenuation) 
measurements. Finally a Tektronix Dual Beam Oscilloscope Type 555 was 
used to display the detected signal and compare it with a monitor signal 
taken directly from the TAB-7 transmitter. Both signals were subject to 
various phase delays (due to coupling into and out of the plasma, ele- 
ments in the detector circuit, etc.), but comparison of the phases of 
the two signals appearing on the oscilloscope gave the relative phase 
delay of the arriving signal. Moving the optical fiber bundle along the 
plasma column until the phase of the arriving signal was the same as be- 
fore allowed a direct determination of the wavelength of the propagating 
signal . 

Signal strength after amplification and filtering was of the order 
of 100 mv (depending on dynode voltages and transmitter power). Signal- 
to-noise ratios varied widely depending primarily on the electronic 
filter settings and on the transmitter power. Individual filter settings 
at each frequency employed were necessary in order to discriminate 
against the noise background in the plasma. The frequencies selected 
for study were chosen in the quieter regions of the plasma noise spec- 
trum. There was a marked variation in the location of these regions 
depending on the static magnetic field strength. A trial and error 
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approach proved to be the most productive method for finding their loca- 
tion. 

A typical oscilloscope photograph is shown in Fig. 12. The upper 
trace shows the detected signal (after amplification and filtering). 

Part of its amplitude is due to an unavoidable inductive pickup condi- 
tion. This part was, of course, subtracted from the raw data before 
analysis. Some plasma noise can be clearly seen within the signal en- 
velope. The s ignal-to-noise ratio here is slightly greater than 1. The 
lower trace is a comparison signal taken from the TAB-7 transmitter. 

This photograph was taken for a 250 kc signal. Deflection sensitivity 
of both traces was 100 mnsr/cm, and sweep speed was 5 ^sec/cm. 

Wavelength and amplitude measurements were taken as functions of 
frequency for various magnetic field strengths. The results provide 
the raw data from \diich the dispersion and attenuation curves presented 
in the next section \jere derived. 

To insure that modes other than theyA= 0 x^re not excited, probe 
measurements were made in an attempt to detect the presence of the 
ylA = 1, 2, 3 modes. (The form of these modes is sketched below.) 

© 

At maximum detection sensitivity no such waves were detected. Since the 
excitation mechanism had azimuthal symmetry, modes other than 0 

should not be present. Actual confirmation that the detected wave was 
actually the U = 0 mode was made by probe measurement. The perturbation 
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Figure 12. Typical Oscilloscope Trace 
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was observed to have a maximum at r = 0 (on the beam axis) and was ob- 
served to decrease as the probe was moved out from the axis. This be- 
havior would only be possible for the ^=0 mode* All other modes would 
exhibit a null at r = 0. This probe measurement also forms the basis for 
the comparison theory for the cutoff distance (that isjZ). These mea- 
surements showed no detectable perturbations beyond a distance of 15 mm 
from the beam axis. 

No standing waves were observed to exist. This was probably due 
to the effect of the burial chamber surrounding the floating anode. In 
the burial chamber ions and electrons are recombined and the resulting 
neutral gas pumped away. Thus the plasma charge density is low, and 
the neutral gas density is high. Both of these conditions would tend 
to preclude the establishment of a standing wave pattern. Reflection 
of the waves from the floating anode could produce a standing wave pat- 
tern, but this effect has never been observed. 17 
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Section 4 



EXPERIMENTAL RESULTS 



Given a fixed static magnetic field strength, a set of data consists 
of (X and ^ as functions of The raw data are reduced in the follow- 

ing way. Assuming cO is real, let ^ = ^ + i^ ( and p both real). 
First the measured wavelengths are reduced to dimensionless variables 
by defining y = . Knowing the wavelength allows direct calculation 

of the damping decrement,^ . If two intensity measurements are made 
one wavelength apart, — II— = e /^l^l ~ where |z, - I " ^ • Thus 

h - ^ " 

— =— = e ' ^ > from which ^ can be calculated directly. The damping 



decrement is then reduced to a dimensionless variable by defining 

. The frequency is reduced to a dimensionless variable by de- 
fining X = . For a given magnetic field strength the entire set of 

data (that is y and z vs. x) are read into a direct search computation 
program. The 5 adjustable parameters W, , C,^, and S are reduced to 
dimensionless variables by defining 



7 



_a 



O' = — 






€ = 



The computation program simultaneously adjusts these parameters to 
find the dispersion function yielding the least square error when com- 
pared to the data. In addition to performing this calculation, graph- 
ical displays are also presented for y vs. x and z vs . x. The plotted 
squares represent the data points, and the continuous curve represents 
the dispersion relation using the 5 parameters just computed. Results 
are presented from 800 to 2200 gauss in the following curves, the para- 
meter being BoCweb/m^). 
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The following 4 curves present the behavior of (the dimension- 
less drift velocity) and (the dimensionless ion thermal speed) plus 
the non- normalized parameters W and C as functions of the magnetic field 
strength. Independent measurement of W was made at 1000 gauss by study- 
ing parallel propagation, and the results of the data reduction for that 
set of dispersion data gave a value within 12% of that obtained from the 
anti-parallel analysis. This will be discussed more thoroughly in Sec- 
tion 5* However, no theories are presented to explain the behavior of 
either parameter. A program was prepared to fit a least squares straight 
line to the deduced parameters as a first approximation. These lines 
are the continuous curves, and the squares represent the fitted para- 
meters for each magnetic field strength. The standard deviations of 
both slope and intercept for these plots appear in Table 1. The units 
of each of the non-normal ized quantities appear on the individual curves. 



Table 1 



Quant ity 




7 '^ 


< 7 ~= 

c 


W(m/s) 


energy 


std. dev. of 


slope 


.00168 


.0741 


12.2 


.120 


std. dev. of 


intercept 


.000264 


.0116 


1.9 


.0188 


slope 




.00109 


.415 


7.92 


.727 


intercept 




-.0124 


.0614 


-90.0 


.00635 
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The next set of 6 curves represents the behavior of ^ (the dimen- 
sionless effective ion collision frequency for momentum loss), £ (the 
dimensionless radial cutoff distance for the perturbations), and € (the 
dimensionless e-folding parameter for the unperturbed radial plasma den- 
sity distribution) plus the corresponding non-normal i zed parameters ^ , 

, and S • Again the squares represent the data points deduced from the 
data reduction, and the continuous curves represent the theories for the 
behavior of these parameters. These theories are presented in detail in 
Section 5. 

By definition j 



and 



e 

~n. 
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The following set of curves depects the behavior of the dimension- 



less phase velocity, — ^ , and the dimensionless group velocity, 

c Jl, 

— , for each value of static magnetic field strength. Again the 

c 

parameter is the value of (web/m^) . The squares represent the normal- 



60 

TST 



ized data points and the continuous curves represent the values of 
and taken from the dispersion relation using the values of the 

5 parameters appropriate to the particular value of 

The modulation feature of the TAB-7 transmitter was inoperative, 
so no group velocity measurements were made, and there is no data with 
which to test the theory. For this reason no data points appear on the 
group velocity curves. 
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The last 4 curves represent the general behavior of the dispersion 



vs. X. The 1000 gauss magnetic field setting was typical, and these 
plots represent the dispersion function with the 5 parameters having 
values appropriate to that setting. Since the 1000 gauss data points 
have already been presented in the previous curves, the data are omitted 
from these 4 curves. 

The first 3 curves are for the complete quadrant, but the group 
velocity is depicted only for the regions in which the damping is small, 
that is where ion acoustic waves might be useful for communication pur- 
poses. These are also the regions where y 0. In the regions of 



vague, since it no longer necessarily represents the velocity of propa- 
gation of a signal. 18 



function throughout the first quadrant of A 5 ^ space. Plots are pre- 




anomalous dispersion < 0) the concept of group velocity becomes 
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Section 5 



INTERPRETATION OF EXPERIMENTAL RESULTS 

Before proceeding to a more detailed analysis of the experimental 
results, the basic behavior of the dispersion curves and the attenuation 
curves will be discussed. Consider first Fig. 13; here a non-dimension- 



isolate the effects of the various parameters on the basic dispersion 
curve, consider first the case of W = 0 (zero ion drift velocity); 



plasma column); the dispersion relation, equation (17) reduces to: 



where o(^ is the first root of J^, = 2.405... 

This equation is plotted in Fig. 13 for the 4 distinct cases: 
= 0, ^ = 00 (zero magnetic field, infinite medium) 



^ ^ (non-zero magnetic field, finite medium) 

The cases presented in Fig. 13 are elementary and are presented 
as an introductory clarification for the more detailed analysis which 
fo Hows . 

From Fig. 13 it is evident that the presence or absence of a mag- 
netic field has no effect on the dispersion in an infinite medium. The 
reason is that the perturbational velocities Vg and are always zero 
for this case (see equations 4a and 4b). Only is non-zero, and since 
it is not a function of SL (see equation 4c), there is no dispersion. 



less plot of ^ vs. is presented. For simplicity, and in order to 



P= 0 (no ion-neutral collisions andJa = Oi + i^ is real); C = 0 (cold 
ions) ; and 5=0 (uniform radial density distribution i^ the unperturbed 




SI ^ 0 ^ ^ CO (non-zero magnetic field, infinite medium) 

= 0, k ^ (zero magnetic field, finite medium) 
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a (ARBITRARY UNITS) vs CJ (ARBITRARY UNITS) FOR THE 
CASE OF W«0, I/ = 0, C*0, 8*0. 



The term cutoff will be used to describe the propagation condition 
under which the phase velocity of the wave is infinite (or when the real 
part of the propagation constant, o( , is zero). The term resonance will 
be used to describe the condition under which the phase velocity is zero 
( is infinite) . 

For a finite medim the high frequencies (and hence small wave- 
lengths) can easily propagate down the plasma column. As the wavelength 
increases there is a point at which the plasma can no longer transmit 
the signal. This appears as a cutoff at frequency ^ , and 

means that at frequencies less than ^0 meaning propagation is 

impossible. This is a consequence of the location of the boundary sur- 
face (the cylinder of radius^). 

Finally for the case of a finite medium with a non-zero magnetic 
field, there is a resonance at the ion cyclotron frequency. At this 
frequency all the ion momemtum appears as cyclotron motion. Again 

2 2 X 

there is a cutoff, but this time at ^ r, = ( ^ — )'^. At fre- 



quencies the plasma beam cannot transmit the signals, and prop- 

agation is impossible since 

In the following discussion it will be assumed that the dispersion 
relation is exactly the Bessel function of order zero. This is equiva- 
lent to the assumption that the unperturbed plasma column has a uniform 
radial density profile. This will considerably simplify the discussion, 
since the infinite series for the dispersion relation need not be consid- 
ered. For this case the relation = 2.405 holds. In actual practice 
the perturbations differ from pure Bessel functions by no more than about 
10%, so the approximation is justified for discussion purposes. 

From this point on the assumption will be made that the plasma is 
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finite in the radial direction 7 ^ oo ) and that there is a non-zero 
magnetic field 0). Also all variables will be normalized (in- 

cluding the 5 parameters) to the dimensionless values. These definitions 
are repeated here for clarification. 



1 , - 
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The effect of 7 ^ (ion drift) is easily understood, since it intro- 
duces a Doppler shift due to the moving medium. In this experiment *5^ 
and were anti-parallel, so the phase velocities are slightly reduced 
by the drift. This is illustrated in Fig. 14. 

A set of measurements was made with and ^ parallel at a field 
strength of 1000 gauss. Direct measurement of the difference in phase 
velocities for the parallel and anti-parallel cases (not the use of com- 
puter data reduction) gave an entirely independent determination of the 
value of the ion drift velocity. This value was within 12% of that ob- 
tained by the computer data reduction process, and this agreement lends 
strong confirmation to the values indirectly deduced from the measure- 
ments . 



The effect of ion temperature is also easily understood. Since ion 

acoustic waves are density perturbations, they propagate at a velocity 

proportional to ( )'^. Passing from the case of cold ions 

inertia term 

to warm ions will increase the pressure term by a small amount equal to 
the ion pressure. Thus the phase velocity will be higher for a given 
frequency; that is ^ will be reduced. This effect is illustrated in 
Fig. 15. The solid curve represents the case of cold ions, and the 
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FIG. 15 

y (ARBITRARY UNITS) VS X( ARBITRARY UNITS) FOR THE CASE 



dashed curve represents the case of warm ions. 

The effect of ^ is more involved. If ion-neutral collisions are 
considered, there will be a decrease in the wave phase velocity to a value 
smaller than in the collisionless case. ^1^ In the present investigation, 
however, this effect is of the second order. The more important effect 
of collisions, however, is that a damping mechanism has been introduced. 
This means that there is no longer an undamped resonance at x = 1. Rather 
the resonance is diffuse, and the maximum wavenumber is limited due to 
the damping. Waves can also exist in the previously forbidden region, 

O2 



1 < X < 






, although here they are heavily damped. This effect is 



illustrated in Fig. 16. Although attempts to excite and detect waves in 



the region 1< x 



*ToT 



were made, the heavy damping plus the long wave- 



lengths in this region precluded their detection. 
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y (ARBITRARY UNITS) v$ X (ARBITRARY UNITS) FOR THE CASE?;* 
O-=0, ^9^0, €=0. 



Ion acoustic waves are produced by the Coulomb interaction (a long 
range interaction). Propagation is possible only when the ion-neutral 
collision frequency is small compared to the wave frequency. The few 
collisions present disrupt the long-range ordering, and damping increases 
with increasing collision frequency. Typical ion-neutral collision fre- 
quencies for this investigation were less than 500 sec"^, so the conti- 
tion of heavy damping would be encountered only at frequencies much 
less than those studied (plus the heavily damped waves in the region 

Kx 

Due to the apparent invariance of ^ with magnatic field, it is sus- 
pected that the mechanism producing the drift is an electric field in the 
interior of the plasma, rather than a fluid pressure differential between 
anodes. If the ion drift were due to the latter mechanism, its value 
would depend on the magnetic field strength, since the neutral gas in- 
jection pressure varied with the field strength. Unfortunately previous 
electric field measurements were not available, so a confirmation of the 
mechanism which produces the drift cannot be made. 

For ion-neutral collisions, X . = ^ 



ion 



^NgO-e 



~7 h. 
(V 



where X . is the ion mean free path 

ion 

is the neutral gas density 

o 

Cr^ is the collision cross section for the process 



(V2) 2 



is the mean ion thermal speed 



Knowing and N (which is manually adjusted for each magnetic 

c g 

field strength to give optimum machine operation), a straightforward 
calculation gives ^ =CP*P x 246 x 10^ where P is the neutral gas 

o o 

pressure and ^ = . 
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Thus 



= 41.7 X 10 




In the above calculation it is assumed that there is a 100% excess ion 
momentum loss per collision. Tliis assumption is not strictly true, but 
it is certainly reasonable, since a large fraction of the directed mo- 
mentum is lost in any collision. It is also assumed that the ion per- 
turbation velocities are small compared to the mean ion thermal speed. 
Thus the collision frequencies will not depend on the strength of the 
pert urbat ion. 

The above equation appears as the theory on the ion temperature 
curves in Section 4. The value of CT" used was taken from the data re- 
duction program appropriate to the magnetic field strength. 

Since 1 (the ions are much colder than the electrons), the dis- 

persion equation is quite insensitive to changes in^ . This can be 
seen by examining the data points of the ^ curves in Section 4. It is 
considered by the author that the deduced values of ^ (and hence the 
ion energies) should not be interpreted too literally. In the ^ vs. 
curve the deduced data points are rather widely scattered. This inac- 
curacy is reflected not only in the ^ vs. curves, but also in the 
vs. Bq curves, since But it can be seen, at least qualita- 

tively, the observed damping results primarily from ion-neutral colli- 
sions . 

Variations in B^ affect the propagation characteristics (that is 

0<2 2 

^2 " (-£1^ — ^^9 since JTL is of the same order of magnitude as 



. An analysis by Wbng concludes that the transverse boundary 





conditions are important only if the ion cyclotron radius, ^ ^ 
parable to the plasma beam radius. At 




IS com- 



mm, so this criterion is approximately met. 
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Neglecting S (uniform unperturbed column) Wong has shown that 



1 7 . 

1 - TToT-y z 






^ o T 4'>f2 0<o2 



where z' is the derivative of the plasma dispersion function 



:(^) = exp(-';;j^2)|^ -2j'exp^2^^j 



with respect to its argument. 

In the present investigation, 






Thus Z ’ (-j~) = 2-^ + — , 

r ^ 2 ’ 



4)r 



X o^Ad^ yy 2T 



and the effect of the radial boundary 



is important only if 2--^ = — 2-lli — . Since this is indeed the case, 

^ /2 

it must be true that T = 1— — L , or T = 2 ev. 

2 £ 

This is a lower limit for the ion temperature, since any tempera- 
ture appreciably lower would predict that ^ 2 a function of 

This estimate is considerably higher than the results of the data re- 
duction (.05 ev to .17 ev), and the spectroscopic measurements of Booth 
and Kelly* j^ioj For that reason it must be stated that the ion energy 
cannot be determined accurately by this experimental method. It is 
probable, however, that the ion energies lie somewhere between .05 ev 
and 2 ev. The ion temperature has been found to increase with distance 
from the plasma source in a helium plasma, varying in value from about 



[-] 



1 ev at the first anode to about 8 ev several meters downstream. 

The increase of ion energy with increasing magnetic field is be- 
lieved to be due to the increased electron-ion collision frequency in 
the denser beams resulting from higher field strengths. 

It should be noted that volume recombination and ionization would 
increase the damping by producing hot neutrals and cold ions respectively. 
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Also, since most ion momentum is directed along the z axis, elastic 
scattering of an ion through anything other than a small angle would 
appear as a damping effect. Charge exchange would substitute a hot 
neutral and a cold ion for the original hot ion and cold neutral. Al- 
though the mean free paths for charge exchange and ion-neutral collisions 
are of the same order of magnitude, the actual collision rate for charge 
exchange would be somewhat less than that for elastic collisions, since 
not every ion-neutral collision would result in charge exchange. But 
charge exchange would certainly appear as another damping mechanism. 

All these processes probably contribute to the inaccuracies of the above 
theory. 

It should be noted that the observed damping is much less than that 
observed by Little and Jones in a gaseous discharge where the percentage 
ionization is low and there are many more neutrals present, Also in 

such a discharge collisions between the ions and the container walls 
would cause additional damping. Such an effect is not a factor in the 
present investigation. 

The primary effect of £ is to alter the form of the radial per- 
turbation from^o to a more rapidly decreasing function of r. The effect 
of £ on y and z is a complicated second order process; it will not be 
discussed further here. 

Neglecting the second-order effects , £ ) , the equations 4a, 

4b, and 4c become (for = 0): 



MNVq 



MNVj. 



3: 




(iiO - +J1 



2 



(itO - 

(id^- ?)2 +J2.2 
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^JL\ 



MNV, = 

^ (iAD -p) 



In this form the equations clearly illustrate the damping effect of 
If "P is now set to zero, these equations can be solved to give: 



= f_L_ ojn 

A> P ^ 



V, 



v„ 









Ap _^2 _ ^2 



I V I = 

I Z I 

Since i?p (except near cutoff), it is clear that most of the 

O ^ 

perturbational momentum appears in the motion of the ions in the z direc- 
tion. It is for this reason that the waves are called acoustic waves. 

Theory for the £ , curves is derived from the probe analysis 
mentioned in Section 3. Those measurements (which confirmed the pre- 
sence of only the ^ = 0 mode) showed that no density perturbation was 
detectable beyond 15 mm from the plasma beam axis. This distance was 
also found to be independent of the magnetic field strength up to 2200 
gauss. This was a remarkably strong confirmation of the values of 
deduced from the data reduction. It also demonstrates the validity of 
the basic dispersion relation, since this relation is far more sensitive 
to small changes in X than for any other parameter. 

Fig. 17 is reproduced from the wrk of Gall and Oleson jisj ; it can 
be seen that a value of /^ = 10 mm would be a more reasonable estimate, 
at least from the point of view of electron temperature alone. However 
the value of = 15 mm taken from the probe analysis gave a much better 
fit to the data and was used as the theoretical value, since it was 



95 




96 



DISTANCE (mm) 



independently confirmed. 



Theory for the ^ , § curves was taken from the work of Gall and 

Oleson. Although their work did not directly yield values of S as 

a function of the magnetic field strength, it did indicate that this 
quantity is not a sensitive function of the magnetic field strength in 
the range used in the present investigation (800 to 2200 gauss). Fig. 
18 is reproduced to illustrate this point. From Fig. 18 S is seen to 
be equal to 100 m"^ at = 600 gauss. Assuming S is not a function 
of Bq, € 



.303 

B^ 



Table 2 summarizes the results of the present investigation for the 
low, medium, and high field strengths used. 



Table 2 



quantity 


800 gauss 


1600 gauss 


2200 gauss 


wavelength (m) 


.0458 to 1.19 


.0314 to .341 


.0248 to .145 


i/e attenuation dist. 
(m) 


20 to 500 


133 to 1110 


455 to 2500 


phase velocity/c 


1.25 to 26 


1.25 to 9.0 


1.2 to 4.4 


group velocity/c 


.09 to .62 


.1 to .73 


. 1 to .8 


ion drift velocity 
(ms“^) 


90.8 


89.4 


87. 1 


effective ion collision 
frequency for momentum 
loss (s”^) 


342 


510 


431 


ion energy (ev) 


.05 


.13 


. 16 


radial cutoff 
distance (m) 


.0152 


.0147 


.0152 


radial 1/e profile 
parameter (m“^) 


.0100 


.0105 


.0105 
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DISTANCE (mm) 



Since the values of y and z vs. x cover such wide ranges, a tabula- 



tion of these quantities vs. x is not included at this point. Represen- 
tative values can be read from the curves presented in Section 4. 

Table 3 



Parameter 


Measured Value 


Deduced Value 


Iwl 


95.2 m/s 


85.0 m/s 




421 s~l 


480 s~^ 


T 


none 


.092 ev 


i 


15 mm 


15.5 mm 


s 


100 m~^ 


102 m"^ 
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Section 6 



CONCLUSIONS 

In conclusion it may be stated that for the propagation of ion 
acoustic waves where Landau damping is not present the new dispersion 
relation derived in Section 2 is capable of predicting the propagation 
characteristics of these waves quite accurately. The greatest error 
present in the dispersion and attenuation curves* (that is the greatest 
departure between data and theory) is 6%. It was, however, much more 
difficult to arrive at satisfactory theories for the 5 adjustable para- 
meters. The close agreement in the dispersion and attenuation curves, 
however, is believed to result from the unambiguous boundary condition 
specified and from the use of photomultiplier detection (rather than 
probe detection with its inherent errors). 

Where damping processes other than ion-neutral collisions occur, 
the effective ion collision rate for momentum loss may be redefined to 
include all such processes (such as the collision rate of ions with the 
walls of the container). For this case ^ would be the sum of the in- 
dividual collision rates for all processes. Since negligible momentum 
is carried by the electrons, electron momentum loss mechanisms may be 
ignored without appreciable error. 

The externally accomplished excitation and dispersion measurements 
described have been a valuable diagnostic tool for the investigation of 
the plasma characteristics. It was by this method that the values of 
W, ^ , U~ , ^ and S were deduced from the data. This method could be 
even more valuable where internal measurements (such as probe diagnos- 
tics) are impractical. 

It is recommended that future investigations be made to test the 
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theory at very low frequencies (1 kc to 100 kc) . Due to the frequency 
limitations of the TAB-7 transmitter employed, measurements in this 
lower frequency region were not possible. Future investigations in this 
region would provide an additional test of the theory and additional con- 
firmation of the deduced values of the 5 parameters considered. 

It is also recommended that the noise spectrum of the plasma be 
analyzed. This would provide additional information on the naturally 
occurring interior phanomena. 

Investigation of ion acoustic wave propagation in other rare gases 
(helium and neon for example) would also provide additional checks of 
the theory. 
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